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Abstract 

We discuss various estimates of the magnitude of higher- twist corrections to the Bjorken 
sum rule in polarized deep inelastic scattering. 
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Considerable attention has been paid to measurements of the polarized structure function 
gi in deep inelastic scattering. The unexpected EMC JTJ results for the first moment of gf, 
the structure function of the proton, invoked tremendous effort to gain more knowledge about 
polarized scattering. As a result of these intense discussions it became clear that polarization 
phenomena might offer much better opportunities to test QCD than unpolarized experiments. 
Meanwhile SMC || and SLAC |§ experiments have provided data on g™ and much improved 
ones on g\. Mainly due to this experiments the general interest shifted to related questions most 
notably to the problem of Q 2 dependence. It turned out that an estimate of the magnitude 
of higher twist contributions is urgently needed in order to compare the experiments of SLAC, 
SMC and EMC which cover different Q 2 regions This is especially important for the SLAC 
data which correspond to Q 2 around 2 GeV 2 . 

The first moment of g±(x, Q 2 ) at fixed Q 2 is given by || 

px Sl (z,Q 2 ) = ia<°> + ^ (a<« + 4<i< 2 > + 4/«>) + 00 

The above formula does not include perturbative corrections (QCD). Higher order corrections 
were recently calculated in the leading twist approximation up to 0(a|) for non-singlet quanti- 
ties, the Bjorken sum rule and - up to 0(a|) - for the Ellis- Jaffe sum rule ||. The Ellis- Jaffe 
sum rule has a flavour singlet contribution. 

The reduced matrix elements for the twist-3 and twist-4 components of Eq. (HD are defined 

by: 

-ipS\q(0) h a g&° + l P gG aa ] q(S>)\pS) - traces = 
6 L J 

= 2d {2) \- (2p a p l3 S a + 1pPp a S a - p l3 p a S a - p a p a S 13 - p a p a S p - p a p l3 S a ) 
—traces] , (2) 



and 



ipS\q(0)gG a ^q(Q)\pS) = 2m 2 N f^S a (3) 



respectively, where \pS) represents the nucleon state of momentum p and spin S, S 2 = —m 2 



A'- 



Note that renormalisation leads to the scale dependence of the matrix elements S 2 ^ 1 and f( 2 \ 
The appropriate scale is given by the value of Q 2 , or the virtuality of the photon probe. 

The reduced matrix elements S 2 ^ and can be expressed through the second moments of 
the polarized nucleon structure functions g\(x) and g^x) M: 



Jdxx 2 g 2 {x) = i(d(2)_ fl (2) 
J dx x 2 gi(x) = — 



(4) 



While may be obtained from existing data in a straightforward way, experimental determi- 
nation of d^ requires high-precision measurement of the g 2 {x) structure function and is much 
more complicated. 

On the theoretical side the higher-twist matrix elements can be estimated either with the 
help of available non-perturbative methods, such as QCD sum rules or ultimately the lattice 
QCD calculation, or within relativistic quark models, such as the MIT bag model. 
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In the usual approach nucleon matrix elements of local operators can be extracted from a 
three-point correlation function 



n r (p) = i 2 / d*xe^ / d A y(0\T{v(x)O r (y)rj(0)} |0> 



(5) 



which involves an interpolating current n(x) with a certain overlap A between the state created 
from the vacuum by r](x) and the nucleon state 

(0\rj(x))\pS) = Xu(p, S) exp (ipx) . (6) 

The overlap integral can be determined from an additional two-point correlation function 

U(p) = i f d^xe^iOlT {r)(x)rj(0)} |0) . (7) 



In practical application it is often advantageous to consider the ratio of three- and two-point 
correlation functions such that the A-dependence cancels out. 

For QCD sum rule calculations of nucleon properties the standard choice for r)(x) has been 
for a long time the three-quark current introduced by Ioffe || 



Vl ( x ) = u a (x)C^u b (x) ^d c {x)e 



(8) 



which was used in the first QCD sum rules calculation of and by Balitsky, Braun and 
Kolesnichenko (BBK) 0. 

As it can be seen from (H) and (|3|), the operators defining d^ and explicitly contain the 
gluonic field operator. This gluonic field has to be matched by another gluonic field operator. 
Because the three-quark current r\i (x) does not contain the gluonic operator explicitly the desired 
contribution must be generated through the perturbative emission of an additional gluon. The 
amplitude for this process is proportional to the strong coupling constant g. 

Alternatively, one may consider an interpolating field which in addition to three quark fields 
contains the gluon field explicitly. In that case gluon emission from the three-quark configuration 
has a non-perturbative character. 

A possible construction of an interpolating current with such properties was discussed in 
It was noticed that all the nice features of the current in @ are preserved by making it 
nonlocal, e.g. by shifting the <i-quark to the point x + e (and adding the required gauge factor 
to preserve gauge invariance). Expanding two times in e and averaging over the directions in 
Euclidean space e fl e u — > \<5^ u e 2 one arrives at 



rf G (x) = e abc (u a (x)C lfl u b (x)) l5lll a a p \gG a/3 (x)d(x 



(9) 



Finally, projecting out the isospin | component leads to the following proton interpolating cur- 
rent: 

(10) 



where 
and 



v £ d (x) = e abc (u a (x)C^u b (x)) 757^ \G^(x)d(x) 



Vg(x) = e abc 



(u a (x)C ltl d b (x)) 757^ [G a ?(x)u(x) 



(11) 
(12) 
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Table 1: Numerical values of higher-twists matrix elements from QCD sum rules [BL [XT 



the MIT 



bag model |jT^, |I7|| , lattice-QCD and experiment |3], |14|, P~5|] . We have extracted value for from 
the published data as the difference between deuteron and proton measurements 2d[ 



(2) 
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It is appropriate to stress here that the technique of QCD sum rules does not require the use 
of "the best" current from all the possible ones, it is only necessary that the current is not too 
bad in order that the contribution of interest is not suppressed by some special reason. 

The QCD sum rules calculation with the current ( |TUD is reported in Refs. [TT, 



Despite 

significant technical differences between both calculations, the numerical results agree satisfac- 



torily. Note that the numerical value of the theoretical uncertainty quoted in jy], [12[] refers only 
to stability of the sum rule. The additional uncertainty, of order of 50%, stems from assumed 
hypothesis about factorization of higher dimensional condensates. This uncertainty has been 
already accounted for in the error estimates presented in ||. 

The most promising technique to calculate matrix elements of local operators in QCD is 
provided, of course, by the lattice formulation. First results for matrix element have been 
recently obtained in [|13|]. The yet unsolved operator mixing problems make at present the lattice 
calculation of impossible. 

The first experimental measurement of d^ 2 \ due to E143 group, has been reported by S. Rock 



during this conference [TJ|], see also (T5J. The data have been taken at average Q 2 ~ 5 GeV 2 . 

MIT bag model estimates are to be found in |16|, [17| . Note that due to SU(6) symmetry 
of the wave function the bag model predicts identically zero for neutron matrix elements. The 
latter bag model estimate |T7j for d^ was obtained from the former one [0 by the perturbative 
QCD evolution, starting from a very low scale /ib ag ~ 0.4 GeV. 

In Table 1 we have summarized various theoretical and experimental estimates of proton (p) 
and neutron (n) matrix elements which contribute to higher-twist corrections to BSR. 

Note that the leading higher-twist matrix elements discussed above describe fundamental 
properties of the nucleon. The twist-4 operator eq. (fj) is a measure for the contribution of the 
collective gluonic field to the spin of the nucleon. Writing the dual field strength tensor in its 
components we get 



(pS\ - B° A f A + (j A x E A y\pS) = 2m 2 N f^S° 



(13) 



where the quark-current is denoted as j A = —gq'yH q and B a A and E a A are the colour magnetic 
and colour electric fields. In the rest system of the nucleon an analogous relation holds for the 
twist-3 operator which determines d^ 

( P S\2B aJa + (j A x E A y\pS) = 8m 2 N d^S° . (14) 

Knowledge of d^ and then allows to estimate magnetic and electric field contributions to 
the spin separately. Using results of [|ll|, [l^] we find that both colour electric and colour magnetic 
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Table 2: Different estimates of the numerical values of the coefficients in front of 0(GeV 2 / Q 2 ) 
corrections to Bjorken (/i(p — n)) and Ellis-Jaffe (p(p)) sum rules, due to QCD sum rules |TT], |T2] ], 
MIT bag model [ff6l O and extrapolation QX8|] based on Gerasimov-Drell-Hearn sum rule. 



fields in the rest system of the nucleon contribute at the same order of magnitude to the spin 

(pS\gB a vlu\pS) = -(0.07±0.08)m^S ff 

(pS\gB a dU\pS) = (0.188 ± 0.08)m 2 N S a 

(pS\ [(ii7 it) x gE J a \pS) = (0.09 ± 0.08)m 2 N S a 

{pS\ [(dj d)xgE] a \pS) = (0.21 ± 0.08)m 2 N S a . (15) 

Obviously such a result shows that simple phenomenological models motivated as analogy to 
QED are misleading. In any such model one would expect the colour- magnetic term to dominate. 

In Table 2 we have summarised various predictions for the coefficient in front of 0(GeV 2 /Q 2 ) 
correction in Eq.(|l|) for Bjorken (p(p — n)) and Ellis-Jaffe (fi(p)) sum rules. The last column 
corresponds to the extrapolation of [Tlj| based on Gerasimov-Drell-Hearn [|l^, |2(| sum rule. It is 
interesting to observe that contrary to the bag model results, the sum rules predictions point in 
the direction suggested by the Gerasimov-Drell-Hearn sum rule. 

It can be seen that the QCD sum rules and the newest MIT bag model predictions, although 
different in sign, are consistent as far as the magnitude is concerned. Hence, one can expect that 
higher-twist effects contribute about 2-3 % of the measured value of the Bjorken sum rule at Q 2 
around 5 GeV 2 . On the other hand further experimental efforts to determine the power-like Q 2 
dependence of the Bjorken integral Jq dxg[~ n {x, Q 2 ) will contribute a lot to our understanding 
to non-perturbative aspects of QCD. 

In the above summary we have purposely omitted discussion of very interesting conceptual 
problem of uniqueness of the definition of power-suppressed corrections to deep-inelastic sum 



rules in QCD. The excellent review of this not yet fully understood topic can be found in [£1 
and references given therein. 
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